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ABSTRACT 



Context. Fast-moving knots (FMK) in the Galactic supernova remnant Cassiopeia A consist mainly of metals and allow to study 
element production in supernovae and shock physics in great detail. 

Aims. We work out theoretically and suggest to observe previously unexplored class of spectral lines - metal recombination lines in 
optical and near-infrared bands - emitted by the cold ionized and cooling plasma in the fast-moving knots. 

Methods. By tracing ion radiative and dielectronic recombination, collisional /-redistribution and radiative cascade processes, we 
compute resulting oxygen, silicon and sulphur recombination line emissivities. It allows us to determine the oxygen recombination 
line fluxes, based on the fast-moving knot model of Sutherland and Dopita (1995b), that predicts existence of highly-ionized ions 
from moderate to very low plasma temperatures. 

Results. The calculations predict oxygen ion recombination line fluxes detectable on modern optical telescopes in the wavelength 
range from 0.5 to 3 fim. Line ratios to collisionally-excited lines will allow to probe in detail the process of rapid cloud cooling after 
passage of a shock front, to test high abundances of 4+ , O s+ and 6+ ions at low temperatures and measure them, to test existing 
theoretical models of a FMK and to build more precise ones. 

Key words, atomic processes - supernovae: individual: Cassiopeia A - infrared: ISM 



1. Introduction 

The brightest source on the radio sky - a supernova remnant 
Cassiopeia A - is so close and so young that the best instruments 
are able to observe exceptionally fine details in its rich structure 
in different spectral bands. One of the most interesting phenom- 
ena are dense ejecta blobs of the supernova explosion, observed 
as numerous bright optical "fast-moving knots" (FMKs). These 
knots radiate mostly in [O m] doublet near 5000 A and other for- 
bidden lines of oxygen, sulphur, silicon and argon. 

Early observations using 200-inch telescope on Mount 
Palomar (iBaade & Minkowskilll954t iMinkowski & Alledll954t 
Minkowskil ll957l) demonstrated that the emitting knot plasma 
has temperatures about 2 x 10 4 K and electron densities about 
10 4 cm~ 3 . Very unusua l spectra of the fast-moving knots have 
lead [Shklovskii (1968) to suggesti on, later confirmed by de- 
tailed measurements and analysi s (Peimb ert & van den B ergh 
1197 It IChevalier & Kirshnerl [l978l |l979j), that they have ex- 
tremely high abundances of some of heavy elements. Some of 
the knots consist of up to 90% of oxygen, others contain pre- 
dominantl y heavier elements, such as silicon, sulphur, argon a nd 
calcium (IChevalier & Kirshnerlll979tlHurford & Fesenlll996l) . 

It was also discovered that a typical lifetime of a FMK is 
of the order of 10-30 years. Existing knots disappear, but other 
bright knots appear on the maps of Cassiopeia A (hereafter 
Cas A). 



* Tables 5 to 23 are only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat?J/A+A/ 



Detailed Ch andra observations (Hughes et alj 2000; 
iHwang et al]|2 004) revealed presence of similar knots in X-rays, 
being very bright in spectral lines of hydrogen- and helium-like 
ions of silicon and sulphur and showing traces of heavily 
absorbed oxygen X-ray lines. 

Similar bright optical knots are observed in other oxygen- 
rich supernova remnants (e.g., Puppis A, N132D, etc., 
ISutherland & Dopita ( 1995a)), but the case of Cas A is the most 
promising for further investigations. 

The FMK optical spectra were first treated as the shock wave 
emission by IChevalier & Kirshnerl dl978l Il979t) . Subsequently 
it was understood that used shock models, that assumed solar 
abundances, are inappropriate for the case of the FMK plasma, 
as its high metal abundances result in a different shock structure. 

There have been published several theoretical model s de- 
scribin g the shock emission in pure oxygen plasmas dltohl 
1 1981 alibi: lB orkowsk i & Shulll l990) and oxygen-dominated plas- 
mas (ISutherland & Dopital I1995U hereafter SD95). In these 
models, the reverse shock of the supernova remnant encoun- 
ters a dense knot and decelerates entering its de nse medium 
dZel'dovich & Raizerll967tlMcKee & Cowiell975l) to velocities 
of several hundred km/s. Shortly after the shock wave enters the 
cloud, the material ahead of the shock front is ionized by the 
radiation from the heated gas after the shock front. All the theo- 
retical models assume a shock wave propagation with a constant 
velocity through a constant density cloud. 

The optical emission is formed in two relatively thin lay- 
ers: one immediately following the shock front, where plasma 
rapidly cools from X-ray emitting temperatures of about 5 x 
10 6 K down to below 10 3 K; another at the photoionization front 
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before the shock wave. The relatively cold layers having temper- 
atures below 10 s K contribute most strongly to the line emission, 
as approximate pressure equilibrium results in much higher den- 
sities and emission measures of these regions. 

Important detail is that the recombination time is much 
longer than the plasma cooling time. Therefore recombining and 
cooling plasma simultaneously contains ions in vastly different 
ionization stages at all temperatures down to several hundred 
Kelvins - the plasma is in a non-equilibrium ionization state. It 
would be very important to find a way to confirm experimentally 
this prediction of the computational models. 

It is interesting to note that predictions of various theoret- 
ical models differ in many ways. For example, as we show in 
an accompanying article (Docenko & Sunyaev, to be submitted), 
the model predictions of the fine-structure far-infrared line in- 
tensities differ by several orders of magnitude. We also show 
there that the models of SD95 and of iBorkowski & Shulll (l990) 
are the best of available ones in reproducing the far-infrared line 
emission of the fast-moving knots, although each of them is only 
precise up to a factor of several. 

Nevertheless, these two theoretical models still have 
a lot of differences bet ween them. For example, in the 
Borkowski & Shull (1990) models the 4+ and higher-ionized 
species are essentially absent, whereas in the SD95 model the 
6+ is the dominant ion after the shock at temperatures above 
about 5000 K. Unfortunately, 4+ , 5+ and 6+ ions have no 
collisionally-excited lines in visible and infrared ranges, and 
their ultraviolet and soft X-ray lines are almost or completely 
undetectable due to high interstellar absorption on the way from 
the Cas A. 

From large differences between the theoretical model predic- 
tions it is clear that they have lack of observational constraints 
and more diagnostic information in form of various line ratios is 
needed to pin down the true structure of fast-moving knots. Part 
of such information can be obtained from the far-infrared lines. 

In this paper we argue that still more information may be 
obtained from metal recombination lines (RL) in optical and 
near-infrared spectral ranges. These lines arise in the transi- 
tions between highly-excited levels (principal quantum number 
re » 5 — 10) populated by the processes of dielectronic and ra- 
diative recombination. The RLs are emitted by all ionic species, 
including, for example, 4+ and 5+ that have no other lines 
detectable from the Cas A. 

Unfortunately, the collisional excitation process that is much 
more efficient at exciting lines corresponding to transitions be- 
tween lowest excited states, is negligibly weak at excitation 
of high-re levels. Therefore, emissivities of the RLs are nor- 
mally several orders of magnitude weaker than ones of the 
collisionally-excited lines, but, as we show in this article, inten- 
sities of the recombination lines from the FMKs are reaching 
detectable levels. 

For oxygen-dominated knots, using plasma parameters 
(electron temperature, density and emission measure of the gas 
at different distances from the shock wave) from the SD95 model 
with the 200 km/s shock speed we have computed the spec- 
trum in the optical and near-infrared bands accessible by the 
modern ground-based telescopes. The cloud shock speed in the 
Cas A FMKs of about 200 km/s is determined by SD95 from the 
comparison of observed optical spectra with their model predic- 
tions. Although similar analysis based on the Borkows ki & Shull 
(1990) model would also be valuable, the original article does 
not give enough information to allow us to compute the RL 
fluxes (e.g., the ionization state distributions as functions of tem- 



perature cannot be resolved from the plots in the regions signifi- 
cantly contributing to the recombination line emission). 

The lines under consideration are predicted to be 300 - 500 
times weaker than e.g. the [Ora] line at 5007 A. However, 
their fluxes are well above the sensitivity of the best present- 
day telescopes. These fluxes are only 2-5 times weaker than 
the detection limits of observations done in 1970's and later on 
the 200-inch and smaller telescopes jPeimbert & van den Bergh 
fl97l IChevalier & Kirshnerl [19791 iHurford & Fesenl 11996). 
Observations of these lines do not demand high spectral res- 
olution due to velocity spreads and turbulence in the shocked 
plasma, as well as due to intrinsic line splitting. 

The recombination lines discussed below are radiated by all 
coexisting ionized oxygen species, from + to 6+ . Study of 
these lines will permit to observe in details the process of the 
non-equilibrium cooling and recombination of the pre- and post- 
shock plasmas overabundant with oxygen or sulphur and silicon. 

The discussed spectral lines allow to observe the coldest ion- 
ized pre-shock and post-shock regions in the FMKs, that are not 
possible to detect using the ground-based observations by any 
other means. The fine-structure far-infrared emission lines are 
also able to give such information, but demand an observatory 
to be located outside the atmosphere and cannot reach angular 
resolution sufficient to distinguish individual knots. 

Observations of recombination lines of elements other than 
oxygen will allow to improve our knowledge on many yet un- 
known physical parameters characterizing the Cas A ejecta. 

The metal recombination lines may be a promising tool 
also for studies of other oxygen-rich supernova remnants like 
Puppis A, N132D, G292+ 1 . 8, etc . 

The paper structure is following. In the second Section we 
describe our method of the recombination line flux computa- 
tion. In the Section [3] we apply the method to predict the op- 
tical and near-infrared line fluxes from the optical knots in 
Cassiopeia A and compare them with existing observational con- 
straints. Description of the recombination line substructure al- 
lowing to identify the parent ion is given in the Section[4] In the 
Section|5]we discuss the ways to get information about physical 
parameters of the emitting region from the line ratios. Finally, in 
Section|6]we conclude the article. 

2. Computation of recombination line fluxes 

In the process of radiative recombination, especially at low tem- 
peratures (when kT is at least several times below the ionization 
potential), significant part of electrons recombine onto excited 
states. Dielectronic recombination for majority of heavy ions 
populate excited states even more efficiently than the radiative 
recombination. 

In low-density plasmas such recombination onto excited 
states, characterized by quantum numbers {nl), is followed by 
the electron radiative cascade to the ground state and emission 
of multiple photons in the course of the cascade. In case of re- 
combining highly-charged ions, these photons are emitted in mi- 
crowave, infrared, optical and ultraviolet spectral bands. 

In this paper we describe the lines produced as a result of 
such radiative cascade, when both collisional and induced tran- 
sitions of the type {nl) — > («'/') are unimportant. The only non- 
radiative process influencing electron cascade in our model are 
the collisional (nl) — > (re/') transitions, as they have much higher 
cross sections and significantly affect the level populations at rel- 
atively low re. They also are the main reason of changes of the 
recombination line emissivities with electron density. 
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We show in Appendix [A] that the n-changing transitions 
may be neglected for computations of optical recombination line 
fluxes, while still retaining reasonable accuracy of the results 
(better than 20%). 

2.1. Elementary processes 

In the current work we account for the following processes: 

- Radiative recombination (RR). Its level- specific rates 
<?rr(«'; T e ) were computed as described in Appendix lB.il 

- Dielectronic recombination (DR). Its level-specific 
rates 2y<?DR(y, nl; T e ) were computed as described in 
Appendix lB.2l 

- Radiative transitions. Their rates A„/„<// were computed using 
hydrogenic fo rmulae with radial integral expressions from 
Gordonl (ll929l) . This may be not accurate for non-hydrogenic 
ions at I < 3, but electrons recombined to those low-/ states 
mostly transit to n' ss / and do not contribute to the optical 
RL emissiorfl 

- Collisional /-redistribution. Its rates C n i„i> were com- 
puted jisin^_jud^eji_TOllMon app roximation expr essions 
from lPengellv & Seatonl dl964l) and lSummersI d 1977b . in the 
region of their applicability. Outside it (where maximum cut- 
off parameter r max be comes less than the minimum one r m ; n ), 
we approximated the lVrinceanu & Flannervl d200ll) classical 
expressions introducing cut-off at large impact parameters 
equal to r max . Details of this approximation are of minor im- 
portance, since the cross sections are proportional as the r^ ax , 
being therefore relatively small. Energies of (n/)-levels and 
their differences determining the collisional rates in many 
cases, were computed as described in Appendix lB.4l 

Usage of hydrogenic approximation does not allow us to re- 
liably compute emissivities and wavelengths of the lines corre- 
sponding to transitions involving levels n < 4. Therefore we do 
not provide the results concerning these levels. 

Details of the atomic physics and approximations utilized for 
computations of the line emissivities are given in AppendixlBl 

2.2. Cascade and I -redistribution equations 

The equation describing radiative cascade of a highly-excited 
electron, accounting for the collisional /-redistribution, is 



N nl(Zn'<n Z/'=/±l A nt,n>l> + ~Ll'±l C nl,nl>) = 

= n e n+ (<7R R (nZ) + 2 y <7dr(7, «0) 



(1) 



where n e , n+ and N„i are the number densities of electrons, re- 
combining ions and recombined ions with electron on the level 
(nl). 

Analytic solution of this coupled system of equations in- 
volves inversion of large matrices of sizes of up to several tens 
of thousands and is not feasible. 

We solved the problem the following way. Starting from 
some maximum n = n max , relevant for the problem (defined in 
the AppendixlAl. we neglected its cascade population. Knowing 



1 After all computations have been performed, the authors were in 
formed about a more elegant and modern way of hydroge nic radial in 



tegral computation using associated Laguerre polynomials (Mali k et"aD 
d 19911) : note that there is a typograp hical error in one of their equations, 
which is rectified in the appendix of Heng & McCray (2007)). 
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Fig. 1. Temperature dependence of the Ov 8 or recombination 
line emissivity in the low-density limit. Lower curve corre- 
spond to DR rate computations using quasi-classical autoion- 
ization rates, upper curve - to the ATOM rates (for details see 
Appendix |B.2| i, middle curve does not account for the process 
of line emission without recombination (DCL), described in 
Appendix IB. 6 1 At low temperatures the emissivity increase is 
determined by the radiative recombination, but the peak around 
10 5 K is arising due to the dielectronic recombination. 



the level-specific recombination and /-redistribution rates, we 
computed resulting populations of levels (nl) by numerical so- 
lution of the system of linear equations (Q]) for n = n max . 

Using level-specific radiative rates, we could compute /- 
resolved cascade population of the levels (nl) for n - n max - 1. 
Having computed total population rates, we solved the sys- 
tem ((TJ for this n. In such a way, moving downwards in n, popu- 
lations of all levels were computed. 

2.3. Recombination line emissivities 

Having obtained level populations N„i, we have computed re- 
combination line emissivities s(nl, n'l'), cm 3 /s, defined as 



e(nl, n'l') 



AN(nl -> n'l') 1 
dVdf n e n 



Iji'I' ■ 



(2) 



where dN(nl — > n'/')/dVdf is the number of transitions from 
level (nl) to (n'l') per second in the unit volume. 

As an example, predicted O v 8a lin^| low-density emissiv- 
ity with /- and /'-components summed up is shown on Figure Q] 
as a function of temperature for different approximations. 

Emissivities of recombining 5+ optical and near-infrared 
recombination lines are given in TableQ]as functions of temper- 
ature in the low-density limit. 

Approximate wavelengths^ of the brightest lines of all ions 
as functions of ionization stage are given in Table [2] They can 
be easily estimated using hydrogenic expressions for any recom- 
bination line arising in transition n — > n' between levels with 
n' > 3, / > 2 as 



A(n — > n) 



0.091127 / 1 



Z 2 



1 



— - — nm, 

n L n L * 



(3) 



2 As it is usual, by net we denote spectral line formed by electronic 
transition from level n + 1 to level n, by nf3 - from n + 1 to n, etc. 

3 Here and everywhere below we present the vacuum wavelengths. 
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Table 1. Wavelengths A and emissivities s of several O v optical 
and near-infrared recombination lines in the low-density limit. 







Emissivity 


e, cm J /s, for electron temperature T c 


Line 


A, yum 


1 • 10 3 K 


3 • 10 3 K 


1 • 10 4 K 


3 • 10 4 K 


5a 


0.298 


9.4 • 10-'" 


4.1 • 10- 12 


3.5 • 10- 12 


5.9 • 10-' 2 


6a 


0.495 


7.6 • 10 12 


3.2- 10~ 12 


1.2 • 10 12 


3.0- 10 12 


la 


0.762 


6.3 • 10-' 2 


2.5 • 10 12 


8.8- 10- 13 


1.9 • 10 12 


8a 


1.112 


5.2 • 10-' 2 


2.0 • 10-' 2 


6.7 • 10~ 13 


1.3 • 10- 12 


9a 


1.554 


4.4 • 10~ 12 


1.7 • 10~ 12 


5.2 • 10~ 13 


8.8 • 10 13 


10a 


2.100 


3.7 • 10~ 12 


1.4- 10- 12 


4.1 • 10- 13 


6.0 • 10" 13 


11a 


2.761 


3.2 ■ 10~ 12 


1.1 ■ 10~ 12 


3.3 • 10~ 13 


4.0 ■ 10" 13 


12a 


3.548 


2.7 • 10~ 12 


9.3 • 10" 13 


2.6- 10~ 13 


2.6 ■ 10~ 13 



Table 2. Hydrogenic vacuum wavelengths, fim, of some optical 
and near-infrared recombination a-lines for several ionization 
stages, denoted by the ion spectroscopic symbol after recombi- 
nation. 



Line 


I 


II 


III 


IV 


V 


VI 


4a 


4.0501 


1.0125 


0.4500 


0.2531 


0.1620 


0.1125 


5a 


7.4558 


1.8640 


0.8284 


0.4660 


0.2982 


0.2071 


6a 


12.365 


3.0913 


1.3739 


0.7728 


0.4946 


0.3435 


7a 


19.052 


4.7629 


2.1168 


1.1907 


0.7621 


0.5292 


8a 


27.788 


6.9471 


3.0876 


1.7368 


1.1115 


0.7719 


9a 


38.849 


9.7122 


4.3165 


2.4281 


1.5540 


1.0791 


10a 


52.506 


13.127 


5.8340 


3.2817 


2.1003 


1.4585 


11a 


69.035 


17.259 


7.6705 


4.3147 


2.7614 


1.9176 


12a 


88.707 


22.177 


9.8563 


5.5442 


3.5483 


2.4641 



where Z is the recombining ion charge. 
2.4. Resulting line fluxes 

Fluxes I(nl, n'l'), erg/cm 2 /s, of the lines were finally computed 
by integrating along the line of sight 



/(«/, n'l') = hv 



— f 

4nR 2 J 



s (nl, n'l'; T e (r)) n e {r)n+(r)dr, 



(4) 



where hv is the photon energy, R is the distance from the ob- 
server to the emitting region and S is the emitting region area. 

The integral over distance in the plain-parallel approxima- 
tion can be easily transformed into integral over temperature by 
substitution 

dr = — — dr = u s hock : dr, 

at dl n x n e AM 

where f s hock is the shock front speed, A# is a cooling function 
and n t is the total number density of all ions in plasma. 

The parameters of this equation - cooling function, electron 
and ion densities as functions of temperature - were taken from 
the SD95 200 km/s shock model, described in more details in 
Section[3]below. 

For purpose of qualitative analysis, we introduce below the 
oxygen differential emission measure per logarithmic tempera- 
ture interval 

d£o rion e dr 



= T e 



dT e 



d(lo g r e ) 

where no is the total oxygen ion number density. It shows con- 
tribution of a given logarithmic temperature interval to the to- 
tal emission measure, showing where most of the line emission 



Undisturbed material 
Photoionization front (T < 30 000 K) 

Cold photoionized region (T ~ 100 K) 



± 



Shock front 



X-ray emitting plasma (T > 10 K) 
Cooling region (T falls to < 300 K~ 
Post-shock photoionized region ? 



Fig. 2. The plain-parallel SD95 model schematic structure. The 
drawing is not to scale. Arrow shows the shock wave direc- 
tion. Optical lines are emitted from the photoionization front and 
the post-shock cooling region at temperatures of several tens of 
thousands Kelvin. Recombination lines arise in the cold pho- 
toionized region and the post-shock cooling region. 



originates. Using this notion, we can express the line flux from 
Eq. © as 



I(nl, n'l') = hv 



4nR 2 



j s (nl,n'l';T e )^^r^T e i\l\. 



d log r e n 



Here the fraction n + /«o is abundance of a given ionic species 
and is also dependent on temperature. 

In our computation of the post-shock recombination line 
emission, we artificially stop integrating expression © when 
plasma temperature drops below 1000 K. This results in some, 
possibly significant, underestimate of the cooling region line 
fluxes arising from recombinations of 2+ ions. Though, as we 
show below, the O n recombination lines from the photoionized 
region are expected to be much brighter. 

It should be remembered that, e.g., O v recombination lines 
arise in the radiative cascade in the 4+ ion triggered by the 
recombination of 5+ ion, so the line fluxes are proportional to 
the ionic abundance of 5+ , not 4+ . 



3. Astrophysical application: FMKs in Cassiopeia A 

The SD95 model describes fast-moving knot emission as arising 
in the interaction of the oxygen-dominated dense cloud with the 
external shock wave, entering the cloud and propagating through 
it (see Figure |2). The heated region just after the shock front 
produces high flux of the ionizing radiation that results in ap- 
pearance of two photoionized regions: before and after the shock 
front. 

The plasma after the shock front passage is rapidly cool- 
ing and, at temperatures of (1 - 5) x 10 4 K, emitting 
the fines observable in the v i sible and near-infrared spec- 
tra (IChevalier & Kirshnerl [l978l fl979t fHurford & Fesenl [1991 
iGerardv & Fesenll2001l) . Thickness of this emitting layer is ex- 
tremely small (about 10 12 - 10° cm), but due to the high density 
of the cooling material at T < 10 5 K it is able to produce bright 
emission lines of highly-ionized atomic species, such as [Om]. 
This region also gives rise to the recombination line (RL) emis- 
sion. 
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The SD95 model does not include the photoionized region 
after the sh ock front. U ntil now, even the presence of this region 
is disputed (Itoh 1986) because, if present, it would produce too 
bright optical lines of neutral oxygen. 

Between the photoionization front and the shock wave the 
plasma is predicted to be extremely cold (T e < 10 3 K) and rather 
highly ionized (ionic abundances of 3+ and 2+ are about 60% 
and 30% in the SD95 200 km/s shock model). This cold ionized 
region, invisible in optical and X-ray bands, should be actively 
recombining and emitting bright recombination lines. 

Resulting RL emission computed based on the SD95 
200 km/s shock model is described separately in the following 
subsections for each of these two regions: after and before the 
shock front. 

In this Section we provide results on the Z-summed recombi- 
nation line fluxes, whereas in the Section|4]we discuss in details 
the individual line spectral substructure, that separates in some 
cases the recombination lines into several components, thus di- 
minishing the individual component fluxes up to factors of 2-4. 

3.1. Cooling region after the shock front 

It is useful to make order-of-magnitude estimates of some 
timescales in the post-shock rapidly cooling plasma at tempera- 
tures between several hundred and 3 x 10 4 K, where the recombi- 
nation line emission occurs (see below). As derived from SD95 
model, the plasma cooling time is roughly equal to 

Tcoq1 30^ w2xl(j4| 

2 n,n e A * ' J 

where T e = 10 4 T 4 K and n t = 10 5 n t ,5 cm . It should be com- 
pared with the time needed to reach the collisional ionization 
equilibrium (CIE), that is in our caseis approximately equal to 
the recombination timescale r rec . For example, for ion 2+ and 
O fi+ the recombination time is 

Trec(0 6+ ) = (ne?RR,0^(re)) _1 * 4.7 X 10 5 S T° M 

T rec (0 2+ ) = (n.q^AT^y 1 * 4.9 x 10 6 s r' ) 64 n e -j, 

where n e = 10 5 n e ,5 cm -3 and <7rr q 2 + is the total radiative re com- 
bination rate for this ion, taken from Pequig not et alj d 199 lb . It is 
easy to see that the recombination timescale in this temperature 
interval is always longer than the cooling time for any abundant 
ion. Therefore the plasma is strongly overionized for its temper- 
ature, that results in high emissivities of the recombination lines. 

The quantitative comparison of these timescales in the model 
is given on Figure [3j where ionization and recombination 
timescales for several ions are compared with the cooling time. 
At high temperatures the cooling timescale is longer than the 
ionization timescales and the material has enough time to con- 
verge to the CIE. At T e < 10 6 K the cooling is much faster 
than the recombination down to temperatures of several hundred 
Kelvin and the high ionization state is stays "frozen" until very 
low temperatures are reached. 

Let us now compare the thermalization timescales a nd mean 
free p aths to the quantities determined above. Following lSpitzerl 
d 1956b . the electron or ion self-collision timescale is equal to 

A U2jV2 



where A is the atomic weight (equal to 1/1836 for electrons), Z is 
the particle charge, n = 10 5 «5 crrT 3 is the particle concentration 




10 3 10 4 10 5 10 6 

Temperature, K 



Fig. 3. Oxygen ion O and O recombination and ionization 
timescales as compared to the plasma cooling time using the 
electron density dependence on temperature and the cooling 
function from the SD95 model with 200 km/s cloud shock speed. 



and A c ~ 3.92x 10 7 n^ /2 J 3/2 . In our case of the highly-charged 
plasma having Z » 1 the timescale of the ion-ion collisions is 
even shorter than for the electron-electron collisions. 

Nevertheless, the mean free path even for electrons A e stays 
significantly below any other characteristic scale, including the 
typical temperature change length scale of more than 10 12 cm: 

A e - T c , e ^3kT e /m e w 1 x 10 7 cm T 2 n~\. 

For illustration of the characteristic conditions after the pas- 
sage of the shock front, on Figure [4] we show Li-like O s+ ion 
density as function of plasma temperature for this regiorQ. The 
figure demonstrates incredibly high abundance of this ion in rel- 
atively dense plasma with T < 3 x 10 4 K after the shock front, 
much higher than in the CIE, where this ion exis ts only in a nar- 
row t emperature range T e = (2 - 7) x 10 5 K ( Mazzotta et al. 
[19981) . 

On the same figure we present also the oxygen emission 
measure distribution over temperature d(EM)/d(logT e ) in the 
SD95 model. This distribution together with the line emissivity 
dependence on temperature (e.g., Fig. [TJ allows one to calcu- 
late the relative contributions to the total line flux from different 
temperature intervals (see Section |2~4l above). 

Such line flux distributions for the prominent [Om] line 
at 5007 A and 8a lines of different oxygen ions are shown 
on Figure [5] Note that at low temperatures (below 4000 K), 
the [Om] 5007 A line emission is dominated by contribution 
from the recombinaiton cascades, not the collisional excitation. 
The cascade contribution to the post-shock intensity of this line 
reaches about 30%. 

Because of high 6+ abundances at both low and high tem- 
peratures, the distribution of the Ovi 8a- line emission is much 

4 We made our own computations of the post-shock plasma recom- 
bination and discovered that our oxygen ion distribution over ioniza- 
tion stages is rather similar to the one presented on the lower panel of 
Figure 3 of SD95, but only if the ion spectroscopic symbols on that fig- 
ure are increased by unity. Therefore we assume that the Figure 3 of 
SD95 paper has a misprint and the ion spectroscopic symbols should be 
read, e.g., "O VII" instead of "O VI", "O VI" instead of "O V", etc. 
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Fig. 4. Differential oxygen emission measure per logarithmic 
temperature interval d£b/d(log T e ) in the post-shock cooling re- 
gion of the SD95 FMK model. Ionic abundances n(O s+ )/n(0) 
of the 5+ ion in the FMK model and in the collisional ioniza- 
tion equilibrium in low-density plasmas (Mazz otta et al.| [l998) 
are also shown. Temperature range not taken into account in our 
calculations is shaded. 




Electron temperature, K 

Fig. 5. Differential luminosity distribution per logarithmic tem- 
perature interval according to the SD95 model for several spec- 
tral lines, including prominent 5007 A [O ra] line and predicted 
On - O vi 8a lines. Curves are normalized so that the area en- 
closed under each of them equals one. Temperature range not 
taken into account in our calculations is shaded. 



wider than that of the collisionally-excited optical [O m] line and 
has 90% of its area within temperature range 5xl0 2 -1.5xl0 4 K. 
Note that the dielectronic recombination is contributing less than 
20% to the resulting line intensity for any of oxygen ions in the 
discussed SD95 model. 

The O v and O vi lines are the brightest expected ones. Lines 
of O in and O iv from this region are predicted to be dimmer by 
a factor of about three, and the O n lines - by a factor of about 
ten. 

On Figure [6] we show the predicted oxygen recombination 
line fluxes from the post-shock region of a model FMK. The 



Table 3. Brightest recombination line /-summed emitted and 
reddened fluxes / and F from the post-shock cooling region 
computed following the SD95 model. The optical [Ora] 5007 A 
line flux also given. The last column contains the flux part in 
the strongest unresolved line component / s . The wavelength A 
corresponds to the strongest line component. 



Ion 


Line 


A, yum 


/, erg/cm z /s 


F, erg/cnr/s 


/. 


Ovi 


7a 


0.5293 


3.43x10" 


16 


3.15x10" 


IK 


0.99 


Ovi 


8a 


0.7720 


1.83x10" 


16 


1.07x10" 


17 


0.99 


Ovi 


9a 


1.0792 


1.03x10" 


16 


2.14x10" 


17 


1.31* 


Ovi 


up 


1.3740 


2.19x10" 


17 


7.42x10" 


18 


0.99 


Ovi 


10a 


1.4586 


6.08x10" 


17 


2.24x10" 


17 


1.00 


Ovi 


13/3 


1.7189 


1.49x10" 


17 


6.65x10" 


18 


1.00 


Ovi 


11a 


1.9177 


3.72x10" 


17 


1.83x10" 


17 


1.00 


Ovi 


12a 


2.4642 


2.35x10" 


17 


1.38x10" 


17 


1.00 


Ov 


6a 


0.4943 


4.07x10" 


16 


2.56x10" 


18 


0.62 


Ov 


7a 


0.7618 


1.99x10" 


16 


1.10x10" 


17 


0.81 


Ov 


8a 


1.1114 


1.06x10" 


16 


2.36x10" 


17 


0.89 


Ov 


10/? 


1.1929 


2.91x10" 


17 


7.66x10" 


18 


0.92 


Ov 


9a 


1.5538 


5.92x10" 


17 


2.36x10" 


17 


1.21* 


Ov 


12/? 


1.9784 


1.24x10 


17 


6.26x10 


18 


0.96 


Ov 


10a 


2.1002 


3.47x10" 


17 


1.83x10" 


17 


0.94 


Ov 


11a 


2.7613 


2.12x10" 


17 


1.32x10" 


17 


0.97 


Oiv 


7a 


1.1902 


2.52x10" 


17 


6.62x10" 


18 


0.80 


Oiv 


8a 


1.7363 


1.32x10" 


17 


5.97x10" 


18 


0.74 


Om 


5a 


0.8262 


6.03x10" 


17 


4.57x10" 


18 


0.25 


Om 


6a 


1.3727 


2.65x10" 


17 


8.95x10" 


18 


0.25 


Om 


7a 


2.1149 


1.28x10" 


17 


6.78x10" 


18 


0.76 


[Ora] 




0.5007 


9.32x10" 


14 


6.30x10" 


16 




[Ora]*" 




0.5007 


1.11x10" 


12 


7.52x10" 


15 





* Including 1 1/? line that is overlapping 
** Including the pre-shock photoionized 



with 9a 

region contribution 



Note. Assumed clo ud area is 1 x 10 33 cm (size about Cf!6). Reddening 
was applied using ISchildl JT977h reddening curve and E(B — V) = 
1.5 dHurford & F esen 1996). The brightest knots are a factor of 3-10 
more intense in the [O in] 5007 A line than the values given in the table. 
The / s values were determined from simulated spectra having 200 km/s 
Doppler line width at the half maximum level (FWH M), that is equal 
to the optical line widths observed in the Cas A knots ( I van den Bergh 
1 197 l\) ; see also Section|4] 



brightest recombination lines around 1 /mi are about 300 times 
less intense than the reddened [O m] 5007 A line. Expected flux 
values are given also in Table [3] 

3.2. Cold photoionized region 

The ions in the photoionized region before the shock front are 
rapidly recombining because of their low temperature. This re- 
sults in rather bright recombination lines, with fluxes propor- 
tional to the distance d between the ionization front and the 
shock wave. 

The plain-parallel SD95 model does not determine this dis- 
tance. Assumption of small optical depth of material between 
the shock wave and the ionization front results in restriction of 
d < 10 17 cm for relevant cloud models for preshock ion num- 
ber density of 100 cm" 3 (illustrated on the Figure 12 of SD95). 
Below we take conservatively of = 10 16 cm, keeping in mind that 
the line fluxes scale linearly with it. 

The cloud model also do not determine exact temperature 
in the photoioniz ed region. We take T e - 1 x 10 3 K as a ref- 
erence value, but Dopita et al] i 19841) mention that it may be as 
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Fig. 6. Reddened recombination line /-summed fluxes F from the 
post-shock cooling region. Details are as in note to Table[3] Here 
the lines are represented as "infinitely thin"; their fine structure 
is discussed in Section|4] 




1.5 2 
Wavelength, \im 

Fig. 7. Recombination line /-summed reddened fluxes expected 
from the photoionized region before the shock front computed 
following the SD95 ionization front model. Details are as in note 
to Table[3] Photoionized region thickness is 10 16 cm. 



Table 4. Brightest recombination line /-summed emitted and 
reddened fluxes / and F from the photoionized region before the 
shock front, computed following the SD95 model. Details are as 
in note to Table[3] Photoionized region thickness is 10 16 cm. The 
last column contains the flux part in the strongest unresolved line 
component / s for FWHM of 200 km/s (see Note to Tabled. The 
wavelength A corresponds to the strongest line component. 



Ion 


Line 


A, yum 


/, erg/cm 2 /s 


F, erg/cnr/s 


/. 


Om 


5a 


0.8262 


4.13xl0" 16 


3.26x10" 


17 


0.26 


Om 


6/3 


0.8327 


1.05xl0" 16 


8.48x10" 


8 


0.31 


Om 


7/8 


1.2552 


6.09x10-" 


1.77x10" 


7 


0.73 


Om 


6a 


1.3727 


1.95X10 -16 


6.61x10" 


7 


0.27 


Om 


8/8 


1.7998 


3.69X10" 17 


1.72x10" 


7 


0.86 


Om 


7a 


2.1149 


l.OOxlO -16 


5.32x10" 


7 


0.82 


Om 


9/8 


2.4804 


2.32x10"" 


1.36x10" 


7 


1.32* 


Om 


8a 


3.0871 


5.50x10"" 


3.59x10" 


7 


0.91 


On 


5a 


1.8645 


4.03x10-" 


1.94x10" 


n 


0.24 


On 


6a 


3.0922 


1.90x10"" 


1.24x10" 


7 


0.59 



* Including lOy line that is overlapping with 9j3 



low as 100 K. If the real temperature in the region is less than our 
assumed value, recombination line emissivities are respectively 
higher (see Section IBTTl below for quantitative dependences). 

Predicted oxygen line fluxes for the SD95 model OP200 with 
shock speed of 200 km/s and total preshock ion number den- 
sity 100 cm" 3 are shown on Figure [7] and given in Table |4] It is 
seen that even with our quite conservative assumptions resulting 
line fluxes are still somewhat stronger than from the shock front. 
Indeed, the predicted line brightest components are on the level 
of almost half a percent of the [O in] 5007 A line. 

If they would be just a factor of several stronger, the lines 
would have been detected in spectroscopic observations of the 
fast-moving knots (see discussion on existing observation limits 
in Section l3~4l below). 

Our estimate of the observational limit on the Om 5a 
recombination line flux at 0.83 /mi based on the data 
bv lHurford & Fesenl d 19961) of 0.005 x F(5007 A) may be trans- 
formed into a joint constraint on the thickness of the photoion- 



ized region d, its temperature T e and total ion density n{. 

I n t \ 2 / 10 3 K \ / d \ 
llOOcm" 3 / \ T e j \10 16 cmj < ' 

assuming E(B - V) = 1.5 and 3+ ionic abundance after the 
ionization front of 0.6. 

On this example it is easy to see that detection of the lines 
or even tight upper limits will allow to perform detailed obser- 
vational tests of the SD95 and other FMK models. 

3.3. Separating pre- and post-shock spectral lines 

Observationally, the lines formed in the regions before and af- 
ter the shock wave will be distinguishable because of two ef- 
fects. Firstly, the lines arising in the photoionized region yet dy- 
namically undisturbed by the shock wave should be intrinsically 
very narrow, but the lines arising in the post-shock region are 
expected to be Doppler-broadened by turbulence. 

Secondly, the spectral lines arising in the post-shock region 
should be shifted with respect to the pre-shock lines. The rela- 
tive Doppler shift between these two line families appears due 
to motion of the post-shock gas relative to the initial knot speed 
induced by the shock wave. In the strong-shock limit this veloc- 
ity difference equals 3/4 1/ s hock for specific heat ratio y — 5/3 
(IZerdovich & Raizerl fl967l) . i.e. up to 150 km/s in the SD95 
model, but the relative Doppler shift will also depend on the an- 
gle between the shock front and the direction towards observer. 

Both collisionally-excited and recombination lines should 
display these features, but for the former group such differences 
should be easier to detect because of higher line intensities. 

The observations with higher spectral resolution (A/AA S= 
20000) could also result in detection of several narrow com- 
ponents forming a line. They could arise from the individual 
vortice or velocity component emission, if there is only one or 
a few of them dominating the plasma motions after the shock, 
like observed in laser experiments simulating cloud -shock in- 
teraction in supernova remnants dKlein et al.l 120031) . The re- 
sulting spectral line profile might be quite complicated (e.g., 
Inogamov & Sunvaev (2003)), but such observations will pro- 
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vide data on the post-shock dynamics of the fast-moving knot 
that is presently poorly constrained from observations. 

3.4. Existing observational limits 

As it is seen from Figures [6] and [7] the brightest recombination 
lines are expected to lie in the wavelength range between 0.7 and 
3 fim. Several detailed spectroscopic investigations have already 
been performed in this range and could have found these lines, 
provided that they had low enough detection limits. Below we 
discuss such existing optical and near-infrared limits on the line 
fluxes and compare them with the SD95 model predictions. 

Optical spectra of fast-mo ving knots with wide spectra l cov- 
erage were obtained, e.g ., by IChevalier & Kirshnerl d 19791) and 
Hurf ord & Fesenl £l996). According to the authors, the detec- 
tion limits around 0.75-0.85 fim are about 300 times less than 
the [O in] 5007 A line flux for the brightest observed features. 

Our results give flux ratios of the reddened 5007 A line to the 
brightest components of the O vi 8a, O v la and O m 5a lines of 
about 300-500 (see Figures and [7] and Tables[3]and|4|. These 
values show that the model is on the boundary of consistency 
with the observational results. They also imply that the oxygen 
recombination line detection should be possible, provided that 
the physics in the post-shock and pre-shock regions corresponds 
to the one described by the SD95 model. 

Recently, near-infrared (0.95-2.5 fim) spectra of Cas A fast- 
moving k nots, obtained at th e 2.4 m Hiltner telescope, were pub- 
lished by Ger ard v & Fesenl (|2001). It is more difficult to com- 
pare their detection limits with our predictions, as there is only 
one oxygen Oi line blend present around 1.129 /urn. It arises in 
transitions between excited states of neutral oxygen and is also 
blended with [S i] line. If we assume that the overlapping [S i] 
line emission is negligible, our estimate of the observational de- 
tection limit corresponds to about 1/100 of the reddened optical 
[O m] line flux. This value is a factor of several higher than the 
predicted fluxes of the O vi 9a, O v 8a and O in 6a lines, also 
showing feasibility of the line detection. 



4. Individual line substructure 

In the non-relativistic hydrogenic approximation, level energies 
depend only on the ion charge and the principal quantum num- 
ber n of the level. Though, this approximation is not fully appli- 
cable for the description of the spectral lines corresponding to 
transitions between the levels with n « 10, as other effects start 
playing a role. They arise due to energy level shifts, described in 
Appendix IB. 41 resulting in the "fine" structure appearing in the 
spectral lines. As a result, these effects help to distinguish lines 
of different elements having equal ionization stages and initial 
and final «'s. 

We calculated the line /- and Zf-substructure following the 
method outlined in Appendix IB. 51 (K is the quantum number 
used to characterize additional interaction of the highly-excited 
electron orbital momentum with the total angular momentum of 
other electrons, present in On, Ora, S n and S m lines). The re- 
sults are shown on Figure [8] for several most interesting cases, 
assuming recombination lines having widths around 200 km/s 
(full width at half maximum level, FWHM), of the same order 
as the observed widths of forbidden optical lines (Ivan den B ergh 
119711) . 

It is easily seen that high-/ levels have both smaller quan- 
tum defects and smaller /f-splitting. Complicated line structure 



in case of presence of the A'-splitting results in weaker individual 
line fluxes making them more difficult to detect. 

One should also keep in mind that spectral lines arising 
in transitions between lower-n levels have higher K- and /- 
splittings. Therefore the line components have larger separations 
and the amount of individual components increase even more 
with each of them being respectivly less intense. As a bottom- 
line, the near-infrared lines are more promising for the first de- 
tection than the optical ones, as they have less substructure (com- 
pare panels (b) and (c) on Figure[8]). 

The substructure of the 6a recombination lines of Mg i, Al i 
and Si i have been observed in solar spectra near 12 fim and ex- 
plained using simil ar, although somewhat more elaborate, theo- 
retical descri ption dBrault & Novesll983l:IChang & Novesll983l: 
I Changl IT984I) . We have used these observations as the test cases 
for our line substructure computation code. 

In Tables 5-15, available as an electronic supplement at the 
CDS, we give O n - O vi, S n - S v, Si n and Si m recombination 
line component vacuum wavelengths, emissivities and relative 
intensities in the (nlK) — > (n'l'K') resolution for several temper- 
atures (lg(r e ,K) =3.0, 3.5, 4.0, 4.5 and 5.0) in the low-density 
limit. First seven columns characterize the quantum numbers 
of the transition (total recombining ion electronic angular mo- 
mentum J c and highly-excited electron quantum numbers before 
and after the transition n, I, K, n' , I', K'), Column 8 gives the 
line component wavelength in microns, next columns in pairs 
state line component emissivity in cm 3 /s and intensity ratio of 
this component with respect to the K- and /-summed emissivity 
e(J c ,nlK,n'l'K')/s(n,n'). Only the lines most likely to be de- 
tected are given in these tables, selected by the following param- 
eters - only a, ft and y lines having wavelengths between 0.3 and 
5.0 fim. 

Tables 16-23, also available at the CDS, contain similar in- 
formation on predicted oxygen line component fluxes in the 
SD95 200 km/s shock model. The fluxes from the cooling and 
photoionized regions are given separately in Tables 16-20 (for 
On - Ovi) and 21-23 (for On - Oiv), respectively. First eight 
columns again contain the quantum numbers characterizing the 
transition and the line wavelength, Column 9 lists line compo- 
nent flux in erg/cm 2 /s, Column 10 contains intensity ratio of this 
component with respect to the K- and /-summed emissivity. 

As discussed in Appendix |B.4| the input atomic data are pre- 
cise to about 10% and we cannot expect better precision of the 
resulting line wavelength differences from the hydrogenic val- 
ues. 

Sample region from the total resulting model spectrum is 
shown on Figure [9] It shows variety of the spectral shapes as 
well as illustrates the diminishing of the peak intensity due to 
Zf-splitting for O m lines. 

The line substructure change with density and temperature 
because of changes in the relative populations of the /-states. At 
low temperatures the main population mechanism is the radiative 
recombination, which is populating high-/ levels relatively effi- 
cient. When the electron density increases, /-redistribution modi- 
fies high-/ populations for n greater than about 20, but for n < 15 
the changes are generally smaller. 

At higher temperatures, DR populates relatively lower / 
states having higher quantum defects and low probabilities of cr- 
ime emission. Therefore a recombination line is split into many 
components and its emissivity is relatively low. In this case, /- 
redistribution mostly shift the recombined electrons to higher-/ 
states, simplifying the line profile and significantly increasing its 
emissivity. 



D. Docenko and R.A. Sunyaev: Optical and NIR recombination lines from Cas A knots 



9 



(8,7) -> (7,6) 
(8,6) -> (7,5) 
(8,5) -> (7,4) 
(8,4) -> (7,3) 



2.07 



2.08 



2.09 2.1 
Wavelength, |im 



2.11 2.12 



< 

CO 

E 



0.74 



O III 5<x+6|3 
(1/50) Si III 5oi+6P 
S III 5a+6|3 



lii 



'■ M It . 




0.75 0.76 0.77 0.78 0.79 0.8 0.81 0.82 0.83 0.84 
Wavelength, |im 



O III 7a 
(1/20) Si III 7a 
S III 7a 



06 2.07 2.08 2.09 2.1 
Wavelength, |im 




2.11 2.12 



1 


O V7a 

- S V 7a 




(d> - 


0.8 








0.6 








0.4 






I: 


0.2 








n 


J \a , . , i 







0.752 0.754 0.756 0.758 0.76 
Wavelength, |im 



0.762 



0.764 



Fig. 8. Fine structure of recombination lines. Line Doppler width at half maximum of 200 km/s is assumed, (a): Om la line 
(8, — > (7, / — 1) component structure for / > 4. Sum of all components is set to be equal for each /. Emissivities of the (8, /) — > 
(1,1 + 1) components are much weaker and structure of resulting lines is not shown. (b),(c): Spectral emissivity ds/dA, cm 3 /s/A of 
the simulated O m, Si m and S m 5a + 6/3 and la line structure for low-density case and temperature r e — 2 x 10 4 K. Note that Si m 
lines have been scaled down to fit in the same scale, (d): Spectral emissivity for the simulated O v and S v la line structure in the 
same conditions. Absence of Zf-splitting results in much simpler profiles. 



On Figure [TO] we show the change of the Om la line fine 
structure with temperature and density, illustrating the described 
effects. 

Another effect arising at low temperatures and electron den- 
sities - lower populations of excited J c core states. For O a, O m, 
S ii and S in lines this results in damping of the /^-splitting, as it 
arises only from recombination of excited ions. All the flux from 
(nlK) — > (n'l'K') components in this case will be redistributed 
into the central («/) — > (n'l') components and only the /-splitting 
will remain. 

5. Plasma diagnostics using recombination lines 

5.1. Temperature diagnostics 

Line emissivity dependence on temperature is shown for several 
bright optical and near-infrared oxygen recombination lines in 
the low-density limit on Figure [TTJ Corresponding line wave- 
lengths are given in Table [2] above. This figure allows to com- 
pute the recombination line fluxes for other models of the multi- 
temperature plasma with emission measure distribution different 
from the discussed SD95 model. 



The curves on Figure [TTJ has two distinct regions. At low 
temperatures the emissivity is determined by the radiative re- 
combination and changes smoothly from ion to ion. At tem- 
peratures higher than several tens of thousands Kelvin the di- 
electronic recombination starts to dominate the recombination 
rates (except for the Ovi lines due to absence at low temper- 
atures of dielectronic recombination channels in He-like ion 
O e+ ) and line emissivities of different ions start being deter- 
mined by the lowest excited states of the recombining ion, as 
described in Appendix IB. 21 The temperature of the strong emis- 
sivity rise is proportional to the typical energy of such excited 
states. Amplitude of the rise is proportional to the lowest excited 
state decay rate. 

For the discussed SD95 200 km/s shock model the high- 
temperature region where the DR dominates is important only 
for O v lines, but it is easy to imagine other ionic abundance dis- 
tributions where it will affect also ions in lower ionization stages. 

Dependence of individual O v recombination or-line ratios on 
temperature is shown on Figure [12] Again, the two distinct re- 
gions are seen at low and high temperatures corresponding to 
RR- and DR-dominated recombination. 
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Fig. 9. Model spectra near 0.8 and 2.1 /urn containing recombi- 
nation lines of O n-O vi computed based on the SD95 200 km/s 
shock model. Contributions of cooling and photoionized regions 
are summed up. It is seen that at longer wavelengths the sub- 
structure is less influencing the line profiles. Note that On ly 
line intensity is multiplied by 10 for illustration purposes. Line 
FWHM of 200 km/s is assumed. 



It is easy to notice that in the intermediate temperatures, the 
lines corresponding to transitions at lowest n's become relatively 
brighter. Reason for such dependence, especially clearly visi- 
ble on the 6a/la line ratio, is the following. As follows from 
Appendix IB. 21 at low temperatures the dielectronic recombina- 
tion rate is damped by the factor exp(-&/kT e ). For the low n- 
levels the doubly excited state energy & may be significantly 
lower than the core excitation energy and at low temperatures 
this makes a difference and the dielectronic recombination pop- 
ulates mostly low-n levels, increasing emissivities of low-n re- 
combination lines. 

The described effect is most pronounced for recombina- 
tion lines of O v, Si in and S v, although in all cases the most 
temperature-sensitive line - the a-line formed by the transition 
from the lowest level, onto which DR is possible - is situated 
around 2000-3000 A. Though, as can be seen from Figure Q~2] 
also the next line situated for these ions around 4500-5000 A, 
can also be a useful diagnostical tool, being relatively much 
brighter in a narrow temperature interval. 
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Fig. 10. Variation of the fine structure of the O m la line at 
2.1 fim with temperature and density. Line FWHM of 200 km/s 
is assumed. Note that the low-density curve for T e = 10 s K is 
multiplied by 20 to be visible on this scale. 




10" 

Temperature, K 

Fig. 11. Low-density /-summed emissivities of oxygen ion opti- 
cal and near-infrared recombination lines. Different colors repre- 
sent different ions, plotted line thickness increase with the spec- 
tral line wavelength. 



5.2. Recombination lines as density diagnostics 

Recombination line flux ratios allow in principle to determine 
also density of the emitting region, as their dependences on den- 
sity and temperature differently affect the relative line fluxes. 

As an example, on FigureQ~3]we plot the O v la line emissiv- 
ity as function of the electron density for different temperatures 
between 10 3 and 10 5 K. For better representation, emissivity ra- 
tios to their low-density values are shown. Clearly, the density 
dependences are weak and can probably be neglected in the first 
stage of the qualitative analysis, especially at low temperatures. 

Two distinct regions are seen once again. At low tempera- 
tures where only radiative recombination determines the level 
populations, emissivities decrease with density, but at higher 
temperatures they increase. The difference is explained by dif- 
ferent initial populations of the levels. 
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Fig. 12. Low-density /-summed O v recombination line emissiv- Fig. 14. Increase of O v a-line emissivities with electron density 
ity ratios as functions of electron temperature. at temperature T e = 3 x 10 4 K. 




Log (n e , cm ) 

Fig. 13. Dependence of O v la line emissivity on electron den- 
sity for different electron temperatures T e . Values are normalized 
to the low-density emissivity. 



Dielectronic recombination predominantly populates states 
with / 5; 10 (see, e.g., Appendix IB.2I >. Collisions in this case 
mostly transfer recombined electrons to higher Z's, increasing the 
probability of An = 1 transitions and, therefore, a-line emissivi- 
ties. 

Radiative recombination, especially at low temperatures, 
populates high-/ states much more efficiently. In this case, pop- 
ulations of low-/ levels relatively increase as a result of the /- 
changing collisions, and recombination line emissivities some- 
what decrease. 

On Figure [14] we show the O v recombination line emissiv- 
ities relative to their low-density values as functions of elec- 
tron density at temperature 30 000 K. Lines arising in transitions 
between higher levels increase relatively more, but in absolute 
terms the increase in emissivity is approximately constant, de- 
termined by the /-redistribution on the high-n levels. 



5.3. Recombination lines of other elements 

Our previous analysis mainly concerned the oxygen recombi- 
nation line emission, as the SD95 model allows making quan- 
titative predictions about their fluxes and line ratios to the ob- 
served collisionally-excited lines. For other elements, no model 
results for ionic abundance distribution are available. Thus to 
predict their line fluxes it is necessary to solve a separate prob- 
lem of non-equilibrium plasma cooling and recombination after 
the shock front for dense clouds of different compositions. Such 
analysis is outside the scope of this paper. 

Nevertheless, it seems valuable to provide data on the re- 
combination line emissivities as functions of temperature for the 
most typical elements. Then, from measured line ratios, it will be 
possible to reconstruct characteristic conditions in the emitting 
regions. As two most typical examples for elements composing 
FMKs in Cas A, we concentrate on the recombination lines of 
silicon and sulphur. Approximate wavelengths of their recombi- 
nation lines are given in Table [2] above. 

Figure [TBI presents low-density emissivities of brightest op- 
tical and near-infrared recombination lines of silicon and sul- 
phur ions, expected to exist in the fast-moving knots. Note that 
Si iv and S vi lines are much weaker than of the other ions at 
T e > 10 4 K due to weakness of dielectronic recombination chan- 
nels at these temperatures. Note also that the line emissivities 
in Si ii and Si m ions start to increase sharply already at about 
12000 K. 

Thus these ions are the most sensitive tracers of plasma at 
temperatures between 15 000 and 30000 K, when emissivities 
of O and S ions are yet relatively weak. This is also seen on 
Figure |8jb,c) above, where the emissivities of the lines of Si m 
are a factor of 20-50 higher than that of O in and S m in the same 
conditions. 

Resembling the case of oxygen recombination lines, the den- 
sity dependences are not pronounced. 

5.4. Ratios to collisionally-excited lines 

Comparison of metal recombination line fluxes with the "tradi- 
tional" collisionally-excited lines may also be a useful tool for 
plasma diagnostics. 
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Fig. 15. Low-density /-summed emissivities of several opti- 
cal and near-infrared silicon and sulphur recombination lines. 
Different colors represent different ions, plotted line thickness 
increase with spectral line wavelength. 



As these two types of lines have different origin, their emis- 
sivity dependences on temperature and density are rather dif- 
ferent. For example, collisionally-excited line emissivities expo- 
nentially drop at temperatures below about hv/k, as thermal elec- 
tron energies are not sufficient to excite the ion electronic transi- 
tion. In contrast, recombination line emissivities at low temper- 
atures increase with the temperature decreasing. 

Also a useful property is that, even in the case when the re- 
combination lines are not detected at all, from the limits on the 
ratios to the collisional lines it is possible to put constraints on 
the plasma parameters. 

On the Figure [16] we give ratios to several brightest 
collisionally-excited optical lines. The exponential increase 
of ratios at low temperatures is clearly seen. Collisional 
line emis sivities have been computed usi ng Chianti atomic 
database dDere et alJll997tlLandi et alj|2006l) 

Ratios to the fine-structure far-infrared line emissivities are 
given on Figure [17] As the Chianti database does not allow to 
compute line emissivities down to 100 K, we computed them by 
extrapolating fine-structure transition electronic excitation col- 
lision strengths to low temperatures by a constant, that should 
be reliable to within a factor of two. The collision strength val - 
ues were adopted from calculations bv lLenno"n & Burkd 71991 : 



Fig. 16. Low-density emissivity ratios of recombination lines to 
several brightest optical collisionally-excited lines. 



Blum & Pradlianl (fl992h : iTaval & Guotal (fl999h : iTavail J2000l 
2006T) . 

We note that observations of the far-infrared lines in the 
spectral range from 10 to 100 /zm are impossible from the 
ground. But even from space, observations of these very intense 
lines cannot be performed with angular resolution sufficient for 
resolving individual knots. They can result only in signal inte- 
grated over many individual emitting objects. 

Density dependences of recombination and collisionally- 
excited lines are also different. Forbidden line emissivity starts 
decreasing as l/n e at some critical density, whereas the recombi- 
nation line emissivities may both decrease and increase depend- 
ing on the plasma temperature (see above). 

The best density indicators for relatively low temperatures 
and densities in the fast-moving knots are obviously ratios of 
different optical and near-infrared line to the far-infrared lines 
that have critical densities of the order of 100 - 10 4 cm 3 . In 
the accompanying paper (Docenko & Sunyaev, to be submitted), 
we present such analysis based on existing experimental data 
and compare the results with our predictions based on the SD95 
model. 

6. Conclusions 

We have offered and developed in detail a new method of rapidly 
recombining plasma diagnostics based on measurements of op- 
tical and near-infrared recombination lines of multiply-ionized 
metal atoms. 

As a promising example, we have applied our method to the 
SD95 theoretical model of fast-moving knots in Cassiopeia A 
supernova remnant and computed expected oxygen line fluxes 
from a single FMK of average size and resulting recombination 
line ratios. It turned out that both cold photoionized region be- 
fore the shock front and rapidly cooling region immediately af- 
ter the shock front produce oxygen recombination lines strong 
enough to be observed on modern optical telescopes in the wave- 
length range between 0.5 and 3 yum. 

At shorter wavelengths, two reasons hamper the line obser- 
vations: high absorption in the interstellar medium on the way 
from Cas A and the line splitting into widely separated multiple 
components with consequently lower intensities in each of them. 
The lines na corresponding to transitions between n > 6 levels 
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Fig. 17. Emissivity ratios of recombination lines to several far- 
infrared collisionally-excited lines. Upper and lower panels 
show emissivity ratios in the low-density limit applicable for the 
cold photoionized region (except for the 88.4 fim line) and the 
case of n e = 3 x 10 4 crcT 3 typical for the cooling region after the 
shock front. 



are the most promising, as they are not so strongly split and have 
50-90% of intensity in one or a few narrow components. 

The precision of our RL flux estimates from the Cas A is not 
expected to be better than a factor of several, as inconsistencies 
of similar magnitude are observed between the SD95 model pre- 
dictions and far-infrared oxygen line observations (Docenko & 
Sunyaev, to be submitted). 

Nevertheless, when detected, the recombination lines will 
allow to determine the details of the photoionization and rapid 
cooling processes in the FMKs from the line intensities, inten- 
sity ratios to each other and to forbidden collisionally-excited 
lines and from the recombination line fine structure measure- 
ments. The measurements of the line structure demand higher 
signal-to-noise ratios that can nevertheless be achieved on mod- 
ern telescopes in a reasonable integration time of few hours. 

One very interesting result is that the predicted O v and O vi 
recombination lines, being brightest in the cooling region recom- 
bination line spectra, arise in the temperature range significantly 
below 10 5 K, where there is essentially no O s+ and 6+ ions in 
the collisional ionization equilibrium conditions. 

It would be also very interesting to observe the recombina- 
tion lines from the knots consisting mostly of silicon and sul- 



phur. Existence of high-temperature shock-heated plasma hav- 
ing such chemical composit ion have been pro v ed by X-ray 
spectr al observations (e.g., iHughes et alj d2000h ; iHwang et al.l 
(2004)), but detailed theoretical predictions of cooling and ion- 
ization structure of such clouds have not been developed yet. 

One of the purposes of this article is to attract attention to the 
metal recombination lines radiated by plasma being far out of 
the collisional ionization equilibrium and having highly-charged 
ions at low temperatures. Besides the supernova remnants, such 
plasma producing detectable metal recombination lines may ex- 
ist in planetary nebulae, as well as near quasars and active galac- 
tic nuclei. 

Acknowledgements. We are grateful to L.A. Vainshtein for valuable advices, 
providing results from the ATOM computer code and permission to use it. We 
are also thankful to the anonymous referee for remarks permitting to make the 
paper easier to understand and for a useful reference on computation of hydro- 
genic radiative transition rates. 

Appendix A: Determination of the upper cutoff n max 

Traditionally, the highly-excited level populations are character- 
ized by the so-called departure coefficients b„i, defined as the 
ratio of actual level population N„i to its thermodynamic equilib- 
rium value N*,. With n increasing, at levels n > n\ the collisional 
/-redistribution processes establish equilibrium population over 
/'s, i.e., b n = b„i. The value of n\ depends on the ion, as well as 
on electron temperature and number density. 

For higher n's, rates of collisional transitions involving n 
change start rapidly increasing and from some n > «2 deter- 
mine the highly-excited level populations. Thus the b„ curve, 
itself defined only at n > ti\, has two distinct ranges as a func- 
tion of n: below n — 112 it is determined by recombination and 
radiative processes and above it rapidly tends to unity because 
of the processes relating different n level populations with each 
other and with the continuum: collisional n-redistribution, three- 
body recombination and collisional ionization. As an example 
relevant for our study, Z?„'s of Li-like ion 5+ at temperatures 
1000 and 20000 K and electron density 50000 cm 3 are shown 
on Figure lATl assuming b„ = b n \. 

On this figure, we also show the departure coefficients, com- 
puted accounting only for radiative n — > n' transitions. It is seen 
that such approximation results in population computation errors 
at levels higher than about tii. 

Let us now show that the electrons that recombined to levels 
n > ri2 are contributing relatively weakly to the total recombina- 
tion and it is safe to introduce an upper cutoff n max , neglecting 
population of higher levels. Two cases —b ni >\ and b m < 1 — 
should be discussed separately here. 

In case of b m > 1 not introducing such cutoff and neglecting 
n—*n' collisions may result in some cases in significant overes- 
timate of the recombination rate - up to a factor of several. The 
easiest estimate of the cutoff position is the level where the total 
recombination rate onto n < n max computed without accounting 
for n — > n' transitions equals to the true total recombination rate. 
It is easy to understand that in this case n max m 112, as illustrated 
in the lower panel of Figure [A3] 

In case of b m < 1 the rates obtained neglecting colli- 
sional transitions between high levels will always underestimate 
true recombination rates. Though for temperatures above about 
10 3 K taking upper cutoff value approximately equal to «2 does 
not result in error more than about 20% for relevant ions, as illus- 
trated on the upper panel of the Figure [A2\ Lower cutoff values 
will result in more significant underestimates of the total recom- 
bination rate and the line emissivities. 
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Fig.A.l. Departure coefficients b n of recombining O ions 
computed accounting for all processes (solid lines) and neglect- 
ing collisional n-redistribution, three-body recombination and 
collisional ionization (dashed lines) at electron temperatures 
T e = 1 x 10 3 K and 2 x 10 4 K and density « e = 5 x 10 4 cirT 3 . 
At higher temperature dielectronic recombination dominates and 
significantly increases the level populations. 



In case of lower densities, when n-i is very high, another 
natural limit « max « (100 - 200) arises from general depen- 
dence of the DR rate o n the principal quantum number, see, e.g., 
ISobelmanet"ai] ( ll98ll) . 



Appendix B: Atomic physics for level population 
computations 

B.1. Radiative recombination 

There are two major recombination processes populating high-n 
states in low-density plasmas: radiative and dielectronic recom- 
bination (RR and DR, respectively). In the radiative recombina- 
tion process, a free electron is directly recombining to a bound 
state (nl). A photon is simultaneously emitted, carrying away the 
released energy. 

The radiative recombination level-specific cross sections 
cr RR were computed from hydrogenic photoionization cross sec- 
tions <x PI by applying the detailed balance relation: 



crr(£ -» nl) = 



a 2 Z 2 21 + 1 



v„(v - v„) 



0p T (raZ -» £), 



(B.l) 



where £ is the electron energy prior to recombination, hv n = 
Ry Z 2 /n 2 is the level binding energy, hv — £ + hv„ is the energy 
of emitted photon, Z is the recombining ion charge and a is the 
fine-structure constant. The photoionization cross sections cr Pr 
were c alculated by the computer program of IStorev & Hummerl 
(1199 11) . 

The level-specific radiative recombination rates 
<?rr(«'; T e ) = (<x RR (£ — > nl) ■ v(&)) were then computed 
by numerical integration over the electron Maxwellian energy 
distribution for the temperature T e (here v denote electron 
velocity). 

Using expressions from e.g. Spitzei] d!956l) . it may be shown 
that the electron-electron collision timescale r ee c is always 
shorter than the cooling time r coo i at temperatures below 10 5 K 
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Fig. A.2. Recombination rates q of 5+ for n e = 5 x 
and T e = 1 x 10 3 K (upper panel) and 2 x 10 4 K (lower panel). 
Meaning of solid and dashed lines is as on Figure I A. 1 1 Plot 
shows both rates for a given n (magenta curves) and cumula- 
tive rates for all levels below n (blue curves). Vertical line show 
upper cutoff n max used for this ion in our calculations. 



where the line emission occurs, and electrons have enough time 
to reach the Maxwellian distribution. Roughly, 



n3/2 



T eec = 2.66 s 



0.16 s 



T 3/2 
1 4 

«e,5 ' 



3.92 x 



«e,5lnA c 

where T e - 10 4 74 K, n e = 10 5 n e 5 cirT 3 and Ac 
10 7 nfT] 11 . 

e,5 4 

B.2. Dielectronic recombination 

The DR is a two- step p rocess (for reviews se e, e. 
Beig man et all d 19681): IShorel i 19691) : ISeaton & Storevl ( fl97... 
Sobelman et al. (119811) 1 illustrated in the simplest case by the 
following diagram: 

X z+ (y ) + e f± X (z - l)+ (y, nl) -> X {Z - l)+ (y , nl) + photon. (B.2) 

Free electron e is first resonantly captured by the ion X z+ 
(Z denotes the ion charge) with simultaneous core electron ex- 
citation from the ground state yrj to the state y having excitation 
energy E c . This process is called dielectronic capture. Resulting 
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ion having two excited electrons is unstable with respect to au- 
toionization - the inverse process of the dielectronic capture. 

The recombined ion then either autoionizes, or its excited 
core electron transits back to the ground state yo and emits a 
resonant line photon, thus making the ion stable against the au- 
toionization. 

Therefore the level-specific dielectronic recombination rate 
<?dr(t, nl; T e ) is a product of two factors: dielectronic capture 
rate and doubly-excited state stabilization probability. In the sim- 
ple case described above it is expressed as 

*»(%»/; TV) = (ne 3/2 s ^mst «*(-*)• ^ 

where y and yo are recombining ion ground and excited states, 
g y and go are their statistical weights, A a (y, nl; yo) is the doubly- 
excited state (y, nl) autoionization rate with core electron after 
the autoionization moving to the ground state yo, A c is the core 
state y decay rate to the ground state yo and S is the electron 
energy prior to the recombination, £ = E c - Ry Z 2 /n 2 . 

We have determined the dielectronic recombination rates us- 
ing two methods, which differ by the way they compute autoion- 
ization rates. In the first method, the autoionization rate is ex- 
pressed in terms of the photoionization cross section using the 
dipole approximation for the inter-electronic interaction. Simple 
expression of this cross section and the autoionization rate in 
the n » 1 limit is then obtained using quasi-classical (QC) ap- 
proach dBureveva & Lisitsall2000l) . 

It is known that the dipole a pproximation is not applicable 
for this purpose for / < 3 states (Beigman et al. 1981). Though, 
this is not a serious limitation in the case, when the dielectronic 
capture occurs via excitation of the recombining ion core elec- 
tron without change of its principal quantum number (An c = 0). 
Such transitions are the most important in our case, as relevant 
plasma temperatures are much lower than the An c > transition 
excitation energies. 

The second method is based on the usage of the program 
ATOM dShevelko & Vainshteinlll993l) . It computes autoioniza- 
tion rates by extrapolating below the threshold the electronic ex- 
citation cross section of the recombining ion core transition. This 
extrapolation method is based on the correspondence principle 
and is thus precise in the limit of n » 1. Both in QC approx- 
imation and ATOM approach, the autoionization rates decrease 
as n~ 3 . 

Comparison plots of the DR rates of O s+ forming 4+ with 
recombined electron populating all n-levels and n — 20 level /- 
states are given on Figure |B~T1 It is seen that in the quasiclassical 
approximation, electrons are populating lower I's, but higher n's. 
As can be inferred from Eq. ( IB. 31 ), this corresponds to sharper /- 
dependence of the autoionization rates. Comparison of scaled 
quasiclassical and quantum autoionization rates A a (2p, nl; 2s) x 
n 3 of doubly-excited oxygen ion 4+ is shown on Figure IB. 21 
indeed showing the inferred dependence. 

It is clear from these data that the two models for the 
DR rates will result in significantly different line emissivi- 
ties (see also Figure [TJ. For final results we decided to use 
the A TOM rates, where available, as the theoretical argu- 
ments dShevelko & Vainshtein|[T993l) show that they are more 
precise than the quasiclassical ones. 

For ions 5+ , 4+ , Si 4+ , Si 3+ , Si 2+ , S 6+ , S 5+ and S 4+ we 
used the ATOM rates. For other ions we implemented the quasi- 
classical expressions. 

Comparison of the dielectronic recombination rates with 
other atomic process rates is shown on Figure llOl for the case of 
4+ ion with highly-excited electron on levels n = 50 and 200. 
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Fig. B.l. Comparison of /-resolved (upper panel) and /-summed 
(lower panel) O s+ ion dielectronic recombination rates com- 
puted by the two methods described in the text: ATOM program 
and quasiclassical (QC) approximation. The values are given for 
T e = 2 x 10 4 K, but the curves scale identically with the temper- 
ature as long as T e «: 10 6 K, when An c > core excitations start 
playing a role. The /-resolved rates are given for n = 20. 



It is seen that the dielectronic recombination populates lower I's, 
but is more efficient at high-n levels. In case of sufficiently high 
density the dielectronic recombination to low-/ states is followed 
by rapid redistribution onto much higher Z's resulting in enhance- 
ment of the recombination line intensities. 



B.3. Highly-excited level energies 

In hydrogen atom the energies of (nl) levels are independent on 
/ because of the shape of the Coulomb potential. In atoms with 
several electrons the changes in the potential energy curve intro- 
duce variations of the highly-excited level energies E(nl) with /. 
This dependence is usually parametrized by expression 



E(nl) = Ry 



(n - flnl) 



2' 



where yU„/ denotes quantum defect of the level (nl) and z — Z + 1 
is the ion spectroscopic symbol. The quantum defects [i„i rapidly 
decrease with increasing / and tend to constant with increasing n 
dSeatonlll983h . 
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Fig. B.2. Comparison of scaled autoionization rates of 4+ ion in 
(2p, nl) state autoionizing to 5+ ion in 2s state, A a (2p, nl; 2s) x 
n 3 , 10 13 s _I , in the quasiclassical approximation and computed 



by the program ATOM. 
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Fig. B.3. Elementary process rates for O s+ ion recombining to 
form 4+ with highly-excited electron having n = 50 (upper 
panel) and n = 200 (lower panel) at electron temperature T e = 
2 x 10 4 K. Collisional rates (total for transitions to all I') are 
shown for electron density n e = 5 x 10 4 cm -3 and are negligibly 
low for n = 50 case. 



Table B.l. Atomic parameters used for computations of the line 
substructure. 



Ion 



On 
Om 
Oiv 
Ov 
Ov 



1.20 
0.97 
0.99 
0.89 



1.35 
0.94 
0.74 
0.27 
0.007 



Ion 



Si n 4.33 
Si m 3.79 
Si iv 



11.7 
7.42 
0.20 



Ion 



Sn 
S in 
S iv 
Sv 
S vi 



3.72 
3.22 
2.44 
2.24 



7.36 
4.79 
3.24 
1.28 
0.07 



Note. Values are given in atomic units (a^ for (r~\ and for a c ). 
Column "ion" contains the ion spectroscopic symbol. 



Their values are important in this study for two reasons. 
Firstly, level energies E(nl) and E(n'l') determine transition en- 
ergies E(nl) - E(n'l') and, therefore, spectral line wavelengths. 
Transitions having the same values of n and n' have similar 
wavelengths. Small wavelength differences are in this case de- 
termined by the energy differences within n and «' groups of 
states - by the quantum defects of the levels. 

Secondly, the collisional transition rates C„i A i< are dependent 
on the energy level splitting \E(nl)-E(nl')\, growing as this split- 
ting decreases. Therefore the quantum defects in this case deter- 
mine rates of /-redistribution and recombination line flux depen- 
dence on electron density in plasma. 

In case of I < 2 the quantum defect s for most ions ar e know n 
and have been taken for our study from Ralc henko et al.l (12007^1 
In the next subsection we describe our method of computation of 
level energy shift from hydrogenic values En(nl) = Ry z 2 /n 2 for 
/ > 2 that can be expressed via quantum defects as 

z 2 

AE(nl) = E(nl) - E H (nl) « 2Ry-r W , 

n i 

where we have assumed that p n \ = i.e. that the quantum defect 
does not depend on n and that it is small: fit <k 1. 

The latter assumption is valid in our case of / » 1, when 
the {nl) states are "non-penetrating", i.e., those in which highly- 
excited electron wave function significantly differs from zero 
only outside the atomic core region. This property also allows 
to use hydrogenic expressions for description of such highly- 
excited states. 



B.4. Computation of the level shifts for high-l states 

Follow ing the approach of Watso rTet al.l (Il980h and iDickinsonl 
(Il981h . we account for the non-penetrating (nl) state level shifts 
arising due to two effects: core polarizability and electrostatic 
quadrupole interaction with the core. 

The energy shift induced by the core polarizability is pro- 
portional to the highly-excited electron radial integral (r~. 4 ) and, 
expressing it with hydrogenic formula, equals (Seaton 1983) 



AEJnl) 



a c z 



2(3 - 1(1+ I) /n 2 ) 



a\ n 3 /(/ + \)(2l - \)(2l + \)(2l + 3) 



Ry. 



(B.4) 



Here ao is the Bohr radius and a c is the c ore polarizability, taken 
for releva nt ions from the handbook of iFraga et al.l (Il976h and 
review bv lLundeenl (T2005) and given for reference in Table IB. ll 
The electrostatic quadrupole interaction is absent for core 
states having total electronic angular momentum / c < 1 or total 
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orbital momentum L c < 1. If J c > 1, the highly-excited elec- 
tron interacts with it, splitting the (nl) level into 2J C + 1 compo- 
nents numbered by the quantum number K, where K — J c + I. 
Normally, if the /^-splitting is present, it is lar ger than the level 
shift due to core polarizability dLundeenl2005l) . 

T he quadrupo le shift in this case is given by the expres- 
sion (Chan^mH: 



AEJnl, K) = (-l) l+L < +s < +K (2J c + 1) 



J l t 2 K \(l\\C^\\l)(L c \\C^\\L c ) (r c 2 ) (r-, 3 )2Ry, 



(B.5) 



where S c is the total core spin momentum, {: : :} denote a 6j- 
symbol, (rfj and (f"~; 3 ) are radial integrals in atomic units for 

core and highly-excited electron, respectively, and (/||C (2) |[/)'s 
are red uced matrix elements, given explicitly, e.g., by Sobelman 
(fl979h : 



(l\\C (2) \\D 



1(1+ 1)(2/+ 1) 
(2/ + 3)(2Z- 1)' 



The highly-excited electron radial integra l can be deter- 
mined using hydrogenic expression from, e.g.. lBethe~& Salpeter 
d!957l) . Core radial integrals ^r 2 ^ were computed for all rele- 
vant ions using radia l wave fu nctions produced by the Flexible 
Atomic Code (FAC, |Gu| d2003l) ). Where literature data for the 
(r^J were available dSenlll979h . the differences between them 
and the FAC values did not exceed about 10-15%. The used val- 
ues of the (r 2 "j are given in Table IB. ll 

Expression flB.5t may be approximated in case of P core 
states (L c = 1) and I » 1 by 



AE a 



M) - 0.2 (r c 2 ) {r-J)Ry 



= 0.2 



nH(l+ !)(/+ 1/2) 



Ry. 



(B.6) 

shift of all 
max (n/) and 



where in this case A£ , q max («Z) gives maximum 
the A"-components. Components with both A£ q 
-A£ , q max (n/) shifts are sometimes present. 

Taking into account these effects, the total (nl) state shift for 
estimate of the collisional /-redistribution rates was taken to be 
maximum possible, i.e., 

AE(nl) = A£ q>max («/) + \AE p (nl)\ . 

As larger shifts correspond to smaller collisional cross sections, 
our approach may somewhat underestimate the /-redistribution 
rates. 

We estimate that the resulting uncertainties of level split- 
tings, in case of the ^-splitting present, correspond to the uncer- 
tainties in the electron density determination from the line ratios 
(see Section |5T2| i being within a factor of about two. 



B.5. Computation of the line substructure 

As the ions of different elements have different values of a c and 
{r^j, the recombination line exact wavelengths will depend not 
only on the ionic charge and the principal quantum numbers of 
the transition, but also on the element and the quantum numbers 
/, 1', K and K'. 

From the above expressions determining energy levels shifts 
it is possible to calculate the positions of all line nl — > n'V com- 
ponents nlK — » n'l'K' . Their relative fluxes are determined by 



the so -called line strengths S according to expression dCh ana 
1984) 



S(nlK, n'l'K') oc (2K' + l)(2K + 1) 



K' K 1 
/ V J c 



(B.7) 



if we assume that populations of the (n//T)-states are proportional 
to their statistical w eights equal to (2K + 1). 

It was shown bv lChang| (ll984) that / = 3 levels in Si i are still 
partially penetrating the core. Therefore also in our case the Si 
and S line substructure is described adequately only for / > 3. In 
our paper we give results also for / = 3, but the precision of line 
position predictions in this case is expected to be rather limited. 

In the post-shock plasma of the fast moving knots, the den- 
sity is high enough to establish significant population of lo west 
excited fine-structure sublevels (Smeding & Pottasch 1979) that 
have typical excitation energies corresponding to temperatures 
of several hundred Kelvin. 

In some cases such excited state populations strongly change 
the line substructure. For example, in 3+ ion, the ground state 
2 Pi/2 has J c = 1/2 and quadrupole interaction with it is absent. 
Though, the lowest excited state 2 P3/2 of this ion has J c = 3/2. In 
case the core electron is in the excited state prior to recombina- 
tion, or getting there after the DR process, the highly-excited lev- 
els will have A'-splitting and resulting recombination lines will 
have much richer substructure. 

In our analysis we assume populations of the 7 C sublevels 
proportional to their statistical weights. 



B.6. Line emission without recombination 

We also account for the following process, mostly not resulting 
in a recombination, but contributing, sometimes noticeably, to 
the spectral line emissivities: 



X Z+ (y ) + e 



-> X iZ - l)+ (y, n'l')+ photon 
-> X z+ (j Q ) + e' + photon, 



(B.8) 



i.e., highly-excited electron transition instead of core transition, 
followed by the captured electron autoionization. This process 
may be quite efficient for increasing An = 1 transition emissiv- 
ities, as the radiative transition rates for captured electron may 
be significantly higher than that of the core electron, if the core 
transition rate is low. 

This is exactly the case for An c = transitions, mostly con- 
tributing to the DR at relatively low temperatures. Rate of the de- 
scribed process ("dielectronic capture with line emission", DCL) 
forming a line from transition (nl) — » (n'V) may be easily ex- 
pressed via the rate of the dielectronic recombination as 

A- l n'V 

4dcl(t, nl, n'l'; T e ) = q D R(y, nl; T e ) " '" . (B.9) 

As the autoionization rates are much higher than the radiative 
transition rates within the ion, probability of the DCL is small 
and after the first captured electron transition it autoionizes, if it 
is energetically allowed. 

To produce noticeable effect, highly-excited electron transi- 
tion rates should be at least comparable to the core transition 
rates, as DCL does not produce cascades. Therefore, only transi- 
tions between lower n's contribute significantly to the additional 
emission in lines. In case of O s+ , described process increases the 
DR contribution to the recombination line emissivity by roughly 
5% for the 8a and la lines. From Eq. dB.9t it is clear that this 
fraction does not depend on temperature. 
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B.7. Estimates of the resulting uncertainties 

Given the amount of approximations made to achieve the final 
results - the metal recombination line emissivities and wave- 
lengths - it might be useful to state explicitly resulting uncer- 
tainties of these quantities. 

The least reliable components for calculation of the recombi- 
nation line emissivities are dielectronic recombination rates and 
upper cutoff position. 

For the ions where the ATOM results are available, the low- 
density emissivity and flux uncertainties are expected to be about 
10%. For the ions where the calculations are made using the qua- 
siclassical expressions, the precision is lower and may constitute 
up to 50%. 

The choice of the upper cutoff n max position results in ad- 
ditional emissivity uncertainty of less than 10% for the DR- 
dominated recombinations and emissivity underestimate of less 
than 20% for the RR-dominated case. 

The density dependences of the line emissivities are not 
strong and should be quite reliably modeled with additional un- 
certainties being below 5-10%. 

It should be reiterated that the absolute oxygen line fluxes are 
much less confidently predicted, as their values are based on the 
theoretical FMK model that shows differences when compared 
to the observed optical-to-infrared collisional line ratios up to a 
factor of several (Docenko & Sunyaev, to be submitted). 

The line fine structure uncertainties are coming from the in- 
put atomic data (polarizabilities and radial integrals) limited pre- 
cision. It is estimated to be not larger than 10-20%, thus the line 
fine structure should be quite reliable. 

The recombination line fine structure component relative in- 
tensities are the least reliably predicted as they depend on the re- 
combining ion ground state populations. The results given above 
assume equilibrium populations; this may not be true at low den- 
sities or temperatures. Therefore the observed line structure may 
be significantly different from one expected. From another side, 
observed line structure will allow to put additional constraints 
on the electron density and temperature from the fine structure 
component intensities. 
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